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ABSTRACT 


This  paper  discusses  infinite  server  queues  whose  input 
is  a Phase  Type  Renewal  Process  introduced  by  M.  F.  Neuts 
(Renewal  Processes  of  Phase  Type,  Nav.  Res.  Log.  Quart, 
forthcoming).  The  problems  of  obtaining  the  transient  and 
steady-state  distributions  and  moments  of  the  queue  length 
are  reduced  to  the  solution  of  certain  well-behaved  systems 
of  linear  differential  equations.  Sample  computations  are 
provided  with  as  many  as  ten  phases.  The  paper  contains 
some  useful  explicit  formulas  and  also  discusses  the 
interesting  special  case  where  the  service  time  is  also  of 
phase  type.  The  Phase  Type  Distributions  include  a wide 
variety  of  models  such  as  generalized  Erlang,  hyper- 
exponential (mixtures  of  a finite  number  of  exponentials) 
as  very  special  cases  and  possess  great  versatality  In 
modeling  a number  of  interesting  qualitative  features  such 
as  bimodality. 


CHAPTER  I 


THE  PH/G/®  QUEUE 


1.1  PH-Distributions  and  PH-Renewal  Processes 

The  versatile  class  of  probability  distributions  of 
Phase  Type  ( PH-di stri buti ons ) which  includes  many  well- 
known  models  such  as  generalized  Erlang,  hyper-exponential 
etc.,  as  special  cases  was  introduced  by  M.  F.  Neuts  [5], 
who,  in  a number  of  follow-up  papers  (c.f.  Bibliography  in 
[7]),  also  discussed  their  useful  computational  properties 
in  the  analysis  of  queueing  systems.  In  this  paper  we  shall 
exclusively  deal  with  PH-di stri butions  of  conti nuous  type 
which  are  obtained  as  the  distribution  of  the  time  till 
absorption  in  an  (m+l)-state  continuous  parameter  Markov 
chain  with  infinitesimal  generator 


Q = 


1° 

0 


and  initial  probability  vector  (a,  am+1 ) , 0<am+1<l,  where 
T=(T1-j)  is  a non-singular  mym  matrix  such  that  T^cO  and 
T.  .*0  for  i^j,  and  T°>0  is  an  m-vector  satisfying  Te+T°=0, 

' J 

where  e 1 = ( 1 , . . . , 1 ) . To  avoid  uninteresting  complications, 
we  shall  henceforth  assume  that  am+i=0.  It  is  now  an  easy 
matter  to  see  that  such  a PH-di stri bution  has  c.d.f., 


F(x)=1 -oexp(Tx)e , x>0.  (1.1.1) 

Such  a distribution  is  said  to  be  of  PH-type  with  represent- 
tation  (cx,T). 

For  technical  reasons  we  may  assume,  without  loss  of 
generality  [5],.  that  the  conservative  stable  matrix  Q*=T+T°A°, 
where  A0ssdiag(a^ , . . . ) and  T°®  (T°  , . . . ,T° ) , is  irreducible. 


Example  1.1:  The  generalized  Erlang  distribution  which  is 
the  convolution  of  m independent  exponential  distributions 
with  parameters,  say,  u-j  »,. . . ,um  respectively  has  a 
representation 


T 


(1  ,0,.  . . , 0 ) 

^ 0 0 

0 -U2  V>2  0 

°°0 


Example  1.2:  The  hyper-exponential  distribution  (which  is 
defined  as  a finite  mixture  of  exponential  distributions)  is 
of  PH-type  with  a representation 

aa(aj.  ,...*.<*  ),•  and  T=di-ag(-ui  , . ...  »-vm) , 

where  u^,...,y  are  the  parameters  of  the  exponential 

distributions  forming  the  mixture  and  c»m  are  the 

ratios  governing  the  mixture. 


A. 


In  [6]  Neuts  discussed  renewal  processes  where  the 
underlying  distribution  is  of  PH-type.  The  "PH-Renewal 
Process"  is  obtained  by  considering  the  Markov  Chain  with 
infinitesimal  generator  Q*=T+T°A°  and  initial  probability 
vector  a,  and  a constructive  definition  of  such  a process  is 
that  it  is  obtained  by  resetting  the  original  Markov  Chain  Q 
following  each  absorption  (i.e.,  renewal)  by  performing  a 
multinomial  trial  with  probabilities  a.  For  such  a PH- 
Renewal  Process,  the  times  between  successive  renewals  are 
i.i.d.,  with  c.d.f.  F(*)  given  by  (1.1.1)  which  is  of  Phase 
Type. 

For  later  use  we  let  e.  denote  the  invariant  probability 
vector  of  the  Markov  Chain  Q* , i.e.,  the  unique  (strictly 
positive)  vector  satisfying 

e_Q*=0 , ee=l . 

We  recall  from  [6]  that  the  stationary  version  of  the  PH- 
Renewal  Process  is  obtained  by  starting  the  Markov  Chain  Q* 
with  initial  probability  vector  e_.  Also  it  is  an  easy 
matter  to  verify  that 

0_=-AaT  1 , 

where  A~^=-aT'^e  is  the  mean  of  the  PH-di stribution  F(*) 
given  by  (1.1.1).  In  the  sequel  we  shall  let  0 denote  the 
mxm  matrix  each  of  whose  rows  is  0. 


r ^ 

r * 

1 . 2 The  PH/G/”  Queue 

In  this  paper  we  consider  a GI/G/»  queue  in  which  the 
arrivals  form  a renewal  process  of  Phase  Type.  Such  a model 

% 

will  be  denoted  by  PH/G/~.  We  shall  assume  that  the  under- 
lying PH- d i stri bution  F C * ) has  a representation  (a,T) 
satisfying  the  conditions  set  forth  in  Section  1.1.  The 
c.d.f.  of  the  service  time  will  be  denoted  by  H(*)»  and  it 

oo  ^ 

shall  be  assumed  that  the  mean  service  time  y=/  { 1 -H ( t ) }dt<«. 

0- 

For  the  PH/G/°°  queue  we  show  that  the  problems  of 
obtaining  the  transient  (and  steady-state)  distributions  and 
moments  of  the  queue  length  (i.e.,  the  number  of  customers 
in  the  system)  can  be  reduced  to  the  solution  of  certain 
well-behaved  systems  of  linear  differential  equations.  By 
applying  a classical  result  on  the  asymptotic  behavior  of 
linear  systems  of  differential  equations,  we  prove,  under 
the  assumption  p<«°,  the  existence  of  a steady-st&te  distri- 
bution and  obtain  an  explicit  formula  for  the  steady-state 
mean  queue  length.  By  considering  the  special  case  where 
H(»)  Is  discrete  with  a finite  number  of  atoms,  we  show  that 
In  the  general  case  there  Is  no  hope  of  solving  these 
differential  equations  explicitly  In  a computationally  useful 
form.  Nevertheless,  these  systems  of  differential  equations 
can  be  solved  numerically  with  considerable  ease.  We  discuss 
some  aspects  of  such  computations  and  present  some  sample 
computations  of  the  first  four  moments  of  the  queue  length 
for  models  with  as  many  as  ten  phases. 


Our  computational  examples  illustrate  how  sole  reliance 
on  the  mean  queue  length  may  result  in  serious  errors  in  the 
qualitative  interpretation  of  the  behavior  of  these 


stochastic  systems.  They  also  enhance  the  importance  of 
algorithmic  solutions  using  numerical  methods  to  stochastic 
models  without  which  an  in-depth  analysis  is  seldom  possible. 

Finally,  we  discuss  briefly  the  interesting  subclass  of 
PH/PH/®  queues  - i.e.,  PH/G/®  queues  where  the  service  time 
distribution  is  also  of  phase  type  - and  in  the  course  of  the 
discussion  obtain  an  alternate  proof  of  a theorem  due  to 
Neuts  and  Jain  [8]  governing  the  independence  of  the  number 
of  customers  in  the  different  phases  of  service  in  an  M/PH/® 
queue . 


CHAPTER  II 


THE  BASIC  SYSTEM  OF  DIFFERENTIAL  EQUATIONS 
AND  ITS  ASYMPTOTIC  BEHAVIOR 


2.1  Introduction 

For  the  PH/G/°°  queue  defined  in  Chapter  I we  let  X ( t ) 
denote  the  queue  length  (i.e.,  the  number  of  customers  in  the 
system)  and  J(t)  the  phase  of  the  arrival  process  at  time  t+. 
For  k>0,  i,j=1,...,m  we  let 

G^j(t)=P[X(t)=k,J(t)=j |X(0)=0,J(0)=i],  t>0  (2.1.1) 

and  let  G^t)  denote  the  m*m  matrix  whose  (i,j)-th  entry  is 
G^J(t).  We  also  define  the  generating  function 

oo 

G(z,t)=  E zkG.(t),  I z I <1  , t>0.  (2.1.2) 

k = 0 K 

In  section  2.2  we  derive  the  system  of  linear  differential 
equations  governing  G(z,«)  and  call  it  the  basic  system  of 
differential  equations  for  the  PH/G/°°  queue.  By  considering 
the  case  where  the  service  time  c.d.f.  H(*)  is  concentrated 
on  a finite  number  of  points,  we  show,  in  Section  2.3,  that 
the  explicit  solution  to  the  basic  system  of  differential 
equations  is  of  no  use  whatsoever  for  computational  purposes 
thereby  suggesting  that  In  the  general  case  there  is  no  hope 
of  obtaining  a useful  explicit  formula  for  G(z,t)  or  its 
limit  as  t->-». 


Finally  In  Section  2.4  we  apply  a classical  result  on 
the  asymptotic  behavior  of  linear  systems  of  differential 
equations  to  the  basic  system  of  the  PH/G/°°  queue  and 
establish  (under  our  assumption  that  the  mean  service  time 
u<»)  that  the  process  (X(t),J(t))  Is  ergodic  and  has  a steady 
state  distribution  independent  of  the  initial  conditions. 

Once  again,  we  do  not  have  an  explicit  formula  for  the  limit. 


2.2 


The  Basic  System  of  Differential  Equations 
Let 


Sn(t)={ (u] 


n 

u.^0,  z u . <t) , t>0, 

i-  ! 1- 


n>l  . 


Also  for  (u-j , . . . ,un)cSn(t) , 0<k<n,  let  g^k  ^ (t  ;u^  , . . . ,un  ) 
denote  the  probability  of  obtaining  exactly  k successes  in  a 
sequence  of  n independent  Bernoulllan  trials  with  respective 
probabilities  of  success  l-H(t-u,),  1 -H  ( t-u-j -u2 ) , ..... 
1-H(t-Ui" . . *-un) . Note  that  1 -H(t-u^ - . . . -u^ ) is  the 
probability  that  a customer  who  arrives  at  epoch  u^+...+u^ 
is  still  in  the  system  at  t. 

In  terms  of  the  above  notations  we  now  have 


Lemma  2.2.1: 

For  t>0, 

GQ(t)=exp(Tt)+ 

00 

Z / exp(Tu, ) *ToA0 • exp(Tu  )-T°A°- 

n-1  Sn(t)  1 n 

exp{T(t-u1-...-un)}g^°^(t;u1,...,un)du1...dun 


and  for  k>l 

oo 

G.  ( t ) = S / exp(Tu1)-T°A°-...exp(Tun)-T°A0- 

k n=k  Sn(t)  1 n 

/ l \ 

exp{T(t-u-|-...-un)}g^  /(t;u1,...,un)du1...dun 


Proof:  The  formulae  above  are  got  by  conditioning  on  n,  the 
number  of  arrivals  in  (0,t],and  the  arrival  epochs  u^  , 
u i +u  2 > • • • » u i + . . . +up  , and  applying  the  law  of  total 
probabi 1 i ty  . 

We  are  now  ready  to  state  our  basic  result  as 

Theorem  2.2.2:  For  all  teC(H),  the  set  of  continuity  points 
of  H C * ) » and  |z|^l,  we  have 

G(z,t)=[(T+T°A°)+(z-l ){l-H(t)}T°A°]G(z,t)  (2.2.3) 

with  the  initial  condition 

G ( z , 0 ) = I (2.2.4) 


Proof : (Throughout  this  proof  we  shall  use  prime  to  denote 

derivatives  with  respect  to  t). 

Noting  that 


n i,  /r\  n M 

E z g ' K ( t ; u , , . . ,u  ) = n [H(t-u, - . . -u . ) + z{ 1 -H(t-u, - . . -u  . ) }]  , 
k=0  n 1 n 1=1  11  11 


we  have  from  Lemma  2.2.1  that 

00  . 00 

G (z  , t ) = E zkG.(t)=  E C ( z , t ) , (2.2.5) 

k=0  K n=0  n 

where 


Cg ( z , t ) =exp ( Tt ) 


and  for  n>l 


C n(z,t)= 

/ ( 8 exp(Tu  . ) •T°A° • [H(t-u, - . . -u  . )+z{1 -H(t-u, - . . -a. ) }]V 

S ( t ) V.i  = 1 1 11  1 1 J 


exp[T(t-u-j-.  . - u )]du-j . . .dun 


Now  for  n>2. 


:n(z,t)  = /|^exp(Tu1  )*T°A°*[H(t-u1)  + zn-H(t-u1)}]* 

f f n exp(Tu.)*T°A°-[H(t-u1-..-u.)+ 
Sn_1(t-u1)0  = 2 1 1 1 

z{l-H(t-u1-..-u1)}]^*exp[T(t-u1-..-un)]* 

di^ • • • dun J du^ 

t 

= /exp(Tu-| ) 'T°A°[H(  t-u-j  )+z { 1 -H(  t-u^ ) }]Cn  i (z  ,t-Uj  )du-j 
t 

= /exp[T(t-T)].T0A0[H(T>+zn-H(T)}]Cn  , (z,x)dx 
0 n" 1 


exp(-Tt)*C(z,t)=/exp(-Tx)*T°A°[z+(l-z)H(T)]C_  , (z ,x )dx  . 

n o n" 1 

This,  on  differentiating  with  respect  to  t,  yields 

-TCn(z,t)+Ci;(z,t)  = T0Ao[z+(l-z)H(t)]Cn_1(z,t)  (2.2.6) 

for  all  n>2  and  teC(H).  It  is  easily  verified  that  (2.2.6) 
holds  for  n=l  also.  Now  adding  (2.2.6)  for  n>l  and  using 
(2.2.5)  we  get 

-T[G(z,t)-C0(z,t)]+[G,(z,t)-Cg(z,t)]=ToA°[z+(l-z)H(t)]G(z,t) 


and  since 


Remark : The  system  of  differential  equations  (2.2.3)  with 
the  initial  condition  (2.2.4)  shall  be  called  the  Basic 
System  of  Differential  Equations  for  the  PH/G/“  queue. 

Note  that  the  basic  system  of  differential  equations  for 
the  PH/G/°°  queue  is  of  the  form 

Y'(t)-[A+R(t)jY(t),  (2.2.7) 

where  A is  a constant  matrix  and  R(t)->-0  as  t-*-«.  Regarding 
the  system  (2.2.7)  we  have  the  following 

Proposition  2.2.8: 

a)  If  R ( • ) is  Riemann-i ntegrabl e In  [0,tQ],  then  there 
exists  a unique  solution  of  (2.2.7)  In  [0,tg]  for  any  given 
initial  condition. 

oo 

b)  If  / | | R(r)  | | d-r<»  and  if  all  the  solutions  of 

0 

Y'(t)=AY(t)  are  bounded,  then  all  the  solutions  of  (2.2.7) 
are  bounded  in  [0,°°). 

Proof : For  Part  (a)  we  refer  the  reader  to  Bellman  [1], 
page  165.  Part  (b)  due  to  Di ni -Hukuhara  may  be  found  in 
Cesari  [2],  p.  37,  3.3  (iii). 


We  can  apply  the  above  Proposition  to  the  system 


(2.2.3) , (2.2.4)  to  obtain 

Theorem  2.2.9:  The  basic  system  of  differential  equations 
for  the  PH/G/“  queue  given  by  (2.2.3)  and  (2.2.4)  has  a 
unique  solution  in  [0,®).  Further  this  solution  is  bounded 
in  [0,-). 

Proof : That  the  second  condition  of  Part  (b)  of  Proposition 
2.2.8  is  satisfied  by  the  system  (2.2.3)  follows  from  the 
fact  that  exp[ (T+T°A° ) t]  is  stochastic  for  any  t>0.  The 
rest  of  the  conditions  in  Proposition  2.2.8  follow  easily 
from  our  assumption  that 

oo 

u = /{l  -H(t)  }dt<°°. 

0 

The  Theorem  now  follows  directly  by  specializing  Proposition 
2.2.8  to  the  system  (2.2.3),  (2.2.4). 

Remark : In  the  special  case  of  the  M/G/«  queue,  the  basic 

system  of  differential  equations  reduces  to  a single 
equation 

ft  G (z ,t )= (z-1 ){l-H(t))XG(z,t) 

G(z ,0)  = 1 , 

for,  in  this  case  T=-x,  T°  = X and  a i = 1 . Under  the  assumption 

oo  ^ 

p = /{l-H(t)  }dt<°°,  we  have 
0 

G(z,t)  = exp^-X/{l-H(T)}dT(l-z)J 
showing  that  X(t)  has  a Poisson  distribution  with  parameter 


x/{ 1 — H ( t ) }df.  The  stationary  distribution  of  the  queue 
0 

length  in  this  case  is  Poisson  with  parameter  Ay.  An 
alternate  proof  of  these  results  can  be  found  in  Takacs  [10] 


2.3  A Discrete  Example 

We  now  consider  the  special  case  where  H(*)  is  concen- 
trated on  a finite  number  of  points  0<tQ<t^  < . . . <tk<“>.  To  I 


specific  we  assume 


1 -H ( t ) = / B, 


where 


if  t<tQ 

if  t^  _i <t<t . ; i = l , . . . ,k 
if  tk<t<°° 


1 = 3q>3-|>.  . . > e k > 0 . In  this  case  it  is  easily  verified 


that  the  solution  to  the  basic  system  of  differential 


equations  (2.2.3),  (2.2.4)  is  given  by 


G(z  ,t)  = 


exp[{ (T+T°A° )+(z-l )T°A° }t] 


if  t < 1 0 


exp[{(T+T°A°)  + (z-l )3iT°A°}(t-ti_1 )]  •G(z,ti_1) 


if  t.  i<t<t . , 1 = 1 ....  ,k 


Lexp[(T+T°A° ) (t-tk)] -G(z  ,tk)  if  t>tk 

which  is  in  the  form  of  a product  of  matrix  exponentials. 
Also  noting  that 

exp{  (T+T°A°  ) t}->0  as  t->°°, 

Lim  G(z,t)=0G(z,t.) 
t->°°  K 

0 

= 0 H exp[{(T+T°A°)  + (z-l  )3,T°A°}(t  .-t . ,)] 
j = k J J J - ' 

where  t_^=0  and  the  matrix  product  is  taken  in  the  order 

j = k,  j = k-l  ....  ,j  = 0. 


The  non-commutativity  of  the  matrices  in  the  products 


above  prevents  any  simplification  and  makes  the  above 
"explicit"  formulae  worthless  for  computing  purposes.  The 
above  example  also  shows  that  in  the  general  case  one  has  no 
hope  of  obtaining  any  useful  explicit  formulae. 

2.4  Asymptotic  Nature  of  the  Basic  System  of  Differential 
Equations 

oo 

Under  the  assumption  ]i  = f{  1 -H (t ) }dt<~  we  show  in  this 

0 

section  that  as  t->-<»  the  unique  solution  G(z,t)  to  the  basic 
system  of  differential  equations  for  the  PH/G/«°  queue  tends 
to  a limit  G(z),  a matrix  whose  rows  are  all  Identical.  This 
is  obtained  by  applying  the  following  classical  result  due  to 
Levinson. 

Proposition  2.4.1:  Consider  the  system  of  m linear 
differential  equations 

d^(t) 

-tt—  = [A+R(t)]y_(t).  (2.4.2) 


Assume  that  the  Jordan  Canonical  form  of  the  constant  matrix 


Assume 


Re  v»A+J. <-3<0  for  j = l , . . . ,s 

and  that 

oo 

/ | r.  . (t)  | dt<<°  for  i,j  = l m. 

0 J 

( k 1 

Then  there  exist  m linearly  independent  vectors  y/  '(t); 
k=l,...,m,  each  a solution  of  (2.4.2)  such  that  as  t-*«. 


y.^  (tj-e^  c(k); 

oi 

• 

• 

* 

II 

J* 

e6t  z(k)(t)+o; 

k= t+1  , . . . 

where  Ac^k^  = y^c^k) ; 

k = l i 

Proof : The  lengthy  proof  of  this  Proposition  may  be  found 
in  Levinson  [3],  Theorem  3. 

oo 

Theorem  2.4.3:  Under  the  assumption  y=/{ 1 -H ( t ) }dt<»,  G(z,t), 

0 

the  unique  solution  to  the  basic  system  of  differential 
equations  (2.2.3),  (2.2.4),  converges  as  t-*°°  to  a matrix 
G(z)  all  whose  rows  are  identical.  Further  G(l-)=0. 

Proof : Since  T+T°A°  is  the  infinitesimal  generator  of  an 
irreducible  continuous  time  Markov-Chai n , (T+T°A°)e  = 0. 

Further  0,  as  an  eigenvalue  of  T+T°A°,  has  multiplicity  1. 
Also  any  other  eigenvalue  of  T+T°A°  has  negative  real  part. 


In  short  T+T°A°  satisfies  the  conditions  for  the  matrix  A in 
Proposition  2.4.1.  It  may  easily  be  verified  that  all  other 
conditions  of  Proposition  2.4.1  are  satisfied  by  the  system 
of  differential  equations 


sg(t) 

-ft—  s [(T+T0A°)  + (z-l){l-H(t)}T°Ao]i(t).  (2.4.4) 

By  Proposition  2.4.1,  the  system  (2.4.4)  has  m linearly 
independent  solutions  £i  (•)»•••  »3fl)(  * ) such  that  as  t-*-« 


£1  ( t )+e 

eBt£k(t)-»-0;  k=2  , . . . ,m 


where  0>-e>Re  X.  for  every  eigenvalue  x.  of  T+T°A°  for  which 
J J 

Re  x .<0.  Now  every  column  of  G(z,t)  satisfies  (2.4.4)  whence 
the  i-th  column  of  G(z,t)  converges  as  t-*»  to  c.j(z)e  where 
c..(z)  is  a constant  which  depends  on  z only.  Thus  as  t-*», 
G(z,t)  converges  to  a matrix  G(z)  all  whose  rows  are 
i denti cal . 

Now  we  can  write  (2.2.3),  (2.2.4)  as 
t 

G(z,t)=I+/[(T+T°A°)+(z-l ){l-H(x)}T0A°]G(z,T)dT,  t>0 
0 

and  therefore 

G(z)=I+/[(T+T°A°)+(z-l ){l-H(T)}T0A0]G(z,x)dT. 

0 

Now, 

G(z)=0G(z) 

00 

= 0+(z-l  )0T°A°/{l-H(x)}G(z,x)dx 
0 


(2.4.5) 


Noting  that  G ( 1 ,t ) = exp[ (T+T°A°  )t]  is  stochastic  for  every 
t>0,  we  can  apply  the  Monotone  Convergence  Theorem  to  show 
that  as  z + 1 , 


oo 

/{ 1 -H(t  ) }G  (z  ,t  )d-re  + Pe. 

0 

Thus  from  (2.4.5)  it  follows  that  as  z + 1 , G(z)-*0,  i.e., 
G(l-)=0,  and  the  proof  is  complete. 

Remarks : 

1.  Note  that  the  (i,j)-th  entry  of  G ( 1 - ) is  the  stationary 
probability  that  the  phase  of  the  arrival  process  is  j. 
Clearly  this  must  be  e.,  for,  e_  is  the  Invariant  probability 

J 

vector  of  the  Markov  Chain  Q*  governing  the  phases. 

2.  The  result  above  also  shows  that  the  rows  of  G(z)  define 
a proper  joint  probability  distribution. 

The  following  Theorem  is  essentially  a re-statement  of 
Theorem  2.4.3  in  the  terminology  of  probability  theory. 

oo 

Theorem  2.4.6:  Under  the  assumption  \i  = f{  1 -H  ( t ) } dx<“ , the 

0 

process  { ( X ( t ) , J ( t ) ) : t>0}  is  ergodic  and  has  a stationary 

distribution  independent  of  the  initial  conditions. 

Proof : A trivial  probabilistic  argument  shows  that  even  if 
we  assume  that  X(0)=k,  J ( 0 ) = i for  any  k>0,  l^i<m,  then  we 
would  still  get  the  same  limit  c.(z)  obtained  in  Theorem 

J 

2.4.3  for  the  sum 


E znP[X(t)=n,J(t)=j|X(0)=k,J(0)=i3 

n = 0 


CHAPTER  III 


MOMENTS  OF  THE  QUEUE  LENGTH 


3.1  Introduction 

In  this  Chapter  we  shall  be  concerned  with  obtaining 
the  moments  of  the  queue  length  X(t).  Letting  N(t)  denote 
the  number  of  arrivals  in  (0,t],  we  can  easily  see  that 
0<X ( t ) <N ( t ) a.s.,  and  since  all  the  moments  of  N(t)  exist 
[6],  so  do  the  moments  of  X(t).  Let  ^jk^(t)  denote  the 
k-th  factorial  moment  of  X ( t ) under  the  assumption  X(0)=0, 
J(0)=i.  That  is, 

u{k)(t)=E[X(k)(t) |X(0)=0,J(0)=iJ,  k*l,  l<i*m, 

where  X ^ k ^ ( t ) denotes  the  factorial  product  X ( t ) [X (t ) -1 ] . . . 
[X(t)-k+l].  We  also  let  \/k^(t)  denote  the  m-vector  whose 
i-th  entry  is  ujk^(t),  k>l . 

Applying  the  rules  of  Calculus  "rather  formally"  to  the 
basic  system  of  differential  equations  (2.2.3)  and  evaluating 

g k 

— r-  G(z,t)e  , we  easily  obtain 

azK  z = 1 - 


jr  n(1 )(t)=(T+T0A0)p(1 )(t)+{l-H(t)}T°A°e 


!i(1)(0)  = 0 


(3.1.1) 


and  for  k>2 


(3.1.2) 


^ w(k)(t)=(T+ToA0)ECk)Ct)+k{l-H(t)}ToAoE(k‘1 \t) 

/k) (0)=0. 

The  above  systems  of  differential  equations  can  easily  be 
solved  numerically  to  obtain  the  first  few  moments  of  the 
queue  size.  In  Section  3.4  we  shall  present  some  sample 
computations  of  the  first  four  moments  with  as  many  as  ten 
phases . 

The  rigorous  proof  of  (3.1.1)  and  (3.1.2)  which  involves 
the  application  of  a number  of  well-known  theorems  in  real 
analysis  is  presented  in  the  next  section.  Section  3.3 
contains  a number  of  remarks  on  the  numerical  computation 
of  the  probability  distribution  and  the  moments  of  the  queue 
length.  Finally  in  Section  3.4  we  present  some  sample 
computations  of  the  first  four  moments  of  the  queue  length 
for  models  with  as  many  as  ten  phases. 

3.2  Moments  of  the  Queue  Length 

We  now  state  a few  well-known  results  from  real  analysis 
which  we  will  need  in  the  course  of  proving  (3.1.1)  and 
(3.1.2). 

Lemma  3.2.1:  Let  ( n , 8 , <j> ) be  a measure  space,  and  let 

(X  : ze[a,b]},  where  [a,b]  Is  a finite  closed  interval  of 

1 ! 

R , be  integrable  functions  from  si  to  R . If  exists  on 

dX 

[a,b]  and  |— ^~-|<Y  where  Y is  integrable,  then 


Proof : See  p.  126,  3°,  Loeve  [4], 


Proposition  3.2.2:  Let  be  a measure  space  and  let 

(X2:  ze[a,b]}  be  integrable  functions  from  n to  Rm.  If  for 

dkX  dkX 

all  k>l  , r-  exists  on  [a,b]  and  || rll<Yb  where  Y.  is 

dzK  dzK  " K K 

integrable,  then 

k dkX 

~ — ir  f X7di)  = / r d(j> , k>l  . 

dzK  n z a dzK 

Proof : This  Proposition  follows  readily  by  repeated 
application  of  Lemma  3.2.1  and  mathematical  induction. 

Proposition  3.2.3:  If  f is  Lebesgue  integrable  on  [a,b]  and 
x 

F(x)=/f (x)dr,  xe[a,b], 
a 

then 

a)  F ( • ) is  a continuous  function  of  bounded  variation 
on  [a,b]. 

b)  F'(x)=f(x)  a.e.  on  [a,b]. 

Proof : For  Part  (a)  we  refer  the  reader  to  Lemma  6,  p.  87, 
Royden  [9].  Part  (b)  Is  Theorem  9,  p.  89,  Royden  [9J. 

We  are  now  ready  to  prove 


Theorem  3.2.4:  For  all  te[0,°°), 

£(1)(t)  = (T+T0A0)/1i.(1  )(T)dT  + /{l-H(T)}dxT° 

0 0“ 


(3.2.5) 


( x ) d x (3.2.6) 


and  for  all  k>2, 
,(k) 


plKJ(t)  = (T+TT)/ji(k)(T)d:  + k/n-H(T)}ToA01i 
0 0 


(k-1) 


Proof : From  (2.2.3)  and  (2.2.4)  It  follows  easily  that 
t 

G(z,t)=I+/[(T+T°A°)+(z-l ){l-H(x)}T°A°]G(z,x)dx, 

0 


whence 

y(k)(t)= 


. k t 


9 z k 0 


/[(T  + T°A°)  + (z-l ){l-H(x)}T°A°]G(z  ,x)dxe 


(3.2.7) 


z = l 


For  ze[0 , 1 ] , x e [ 0 , t]  , clearly, 
, k 


0 < 


- vjk)(t) 


< E[N(k)(x)|J(0)=i] 

s E[N{k)(t)|J(0)=i], 

where  N ^ k ^ ( t ) = N ( t ) [N ( t ) - 1 ] . . . [N ( t ) -k+1 ] . The  last  inequality 

in  the  above  chain  of  inequalities  is  got  by  using  the  fact 

that  Nk(*)/^a-s.  Now,  denoting  E(N^  (x)  | J(0)  = ij  by 
( k 1 

vj  ' ( x ) , we  have 
t ( k \ 

0 < I v.  (t)dx<»  for  every  l<i<m,  k>l  . 

0 1 

Also  for  0<z<l , 0<x<t , 


A. 


| |^—y[  {(T+T°A°)+(z-l  ){1-H(t)}T°A0}G(z,t)£J| | 

3k 

= | |[(T+T0A0)  + (z-1){1-H(t)}T0A0J— - G(z,x)e 

3ZK 

,k«l 

+ k{l-H(T)}T0A°  - j-  t-  G(z,r)e|  | 
azK 

<{ | | T + T°A0  | | + | | T°A°  | | } | |v(k)(t)  | |+k|  |T°A°  | | | | v ( k" 1 } ( t ) | | , 

where,  \/k^(t)  is  the  vector  whose  i-th  component  is  v(k^(t). 
Now,  in  Proposition  3.2.2  set  ft=[0,t],  B=Borel  subsets  of 
[0,t],  4>  = Lebesgue  measure, 

Xz(t)=[(T+T°A°)+(z-1 ){l-H(T)}T°A°]G(z,T)e 

Y k ( x ) = ( | | T+T0A°  | | + | | T° A° | | )|  |v(k)(t)|  | +k  | |T°A°|  | | |v(k_1)(t)|  | 
to  obtain,  using  (3.2.7),  that 

i(k)(tW  iLl(T+T»A0)  + (z-l)U-fl(T>n"A"]S(z,T)e|  dt, 

0 Hz-1- 

and  the  theorem  follows. 

Corollary  3.2.8: 

( k 1 

a)  For  all  k > 1 , '(•)  is  continuous  and  of  bounded 

variation  in  any  finite  interval  [0  , t ] . 

b)  We  also  have 

£^^(t)  = (T+T0A0)y^1)(t)  + {l-H(t)}T0A0f^,  a.e. 

Uo)  = o 

and  for  k>2 

u(k)(t)=(T+T0A°)]i(k)(t)  + k{l-H(t)}T0A0v(k"1  J(t)  , a.e. 

^k  ^ (0)  = 0 


Proof:  This  Corollary  is  an  immediate  consequence  of 
Theorem  3.2.4  and  Proposition  3.2.3. 


The  above  discussion  completes  the  proof  of  (3.1.1)  and 
(3.1.2)  which  were  derived  heuri stical ly  at  the  beginning  of 
this  Chapter.  In  a very  easy  way  we  can  now  prove 


Theorem  3.2.9:  Under  the  assumption  y = /{ 1 -H (t ) }dx<» , for 

0 ( k 1 

any  k>l  the  vector  of  factorial  moments  v}  '(t)  converges 
as  t-*°°  to  a (finite)  vector  all  whose  components  are  equal. 
Further  for  any  k > 1 , y^k^(*)  is  bounded  in  [0,°°). 


Proof:  Clearly  it  suffices  to  prove  the  existence  and 

finiteness  of  these  limits.  To  this  end  note  that  (3.2.5) 
implies 

0y  ^ 1 ^ (t  ) = 0T°/{  1 - H (t  ) }dT-nteT°  <°°  as  t-*». 

0 

This  shows  that  y^^(t)  converges  to  a finite  limit  as  t-*-~ 
and  also  (since  y.-|  ( • ) is  continuous  in  [0,°°))  that  ^ ( • ) is 
bounded  in  [0,°°).  Now  as  induction  hypothesis  assume  that 
the  Theorem  is  true  for  k-1 . From  (3.2.6)  we  now  have 

0y(k)(t)  = 0ToAo/{l-H(T)}y(k"1  ^xjdx, 

0 

and  the  above  integral  converges  as  t+“>  to  a finite  limit 
since  y/  '(•)  is  bounded  on  [0,«>).  Now  the  boundedness  of 
y (•)  follows  from  its  continuity.  By  mathematical  induction, 
the  proof  is  complete. 


A 


A 


An  interesting  special  case  of  the  PH/G/«°  queue  is  the 
one  in  which  the  PH-Renewal  process  describing  arrivals  is 


stationary  - i.e.  where  the  initial  phase  is  chosen  according 
to  the  vector  For  this  case  we  have  the  interesting 


Theorem  3.2.10:  Suppose  the  initial  phase  J ( 0 ) is  chosen 
according  to  the  vector  e_,  i.e.  the  PH-Renewal  process  is 
stationary.  Then 

t 

a)  the  mean  system  size  at  t is  given  by  \f{ 1 -H (t  ) }dr . 

0 

co 

b)  if  p = /{ 1 -H ( x ) }dT<» , then  the  mean  of  the  steady 
state  distribution  of  the  system  size  is  \v  where  \ is  the 
mean  inter-arrival  time. 

Proof : Pre-mul  ti plyi ng  (3.2.5)  by  e_, 

eu(1)(t)=eT0/{l-H(T)}dr, 

0 

and  the  result  in  Part  (a)  follows  by  noting  that  e^-AoT"^ 
and  T°=-Te.  Part  (b)  now  easily  follows  from  Part  (a). 

Using  the  above  theorem  and  Theorem  3.2.9  we  immediately 
obtain 

oo 

Corollary  3.2.11:  If  y = /{  1 -H  ( t ) }dx<°° , then 

0 

^ ^ ( t ) -*A y e^  as  t -»■<». 

3.3  Some  Remarks  on  Computational  Methods 

Except  in  the  special  case  of  Poisson  arrivals,  there 
does  not  seem  to  be  any  hope  of  obtaining  in  closed  form  the 
distribution  of  the  system  size  or  the  moments  thereof  for 


A 


A 


the  PH/G/<«  queue.  Nevertheless,  all  the  systems  of 
differential  equations  above  lend  themselves  readily  to 
computations  using  numerical  methods. 


m 


For  the  purpose  of  computing  the  moments  up  to  an  index 
k,  it  appears  best  to  consider  the  km  differential  equations 
given  by  (3.1.1)  and  (3.1.2)  as  forming  a single  system. 
Using  a general  purpose  Inter-active  software  system  called 
DELSIM,  due  to  Professor  D.  E.  Lamb  and  available  at  the 


University  of  Delaware  Computing  Center,  we  computed  the 
time-dependent  solutions  to  the  moment-vectors  up  to  the 
fourth  moment  for  models  involving  as  many  as  ten  phases  by 
applying  the  Fifth-order  Kutta-Merson  method.  Some  of  these 
computations  are  presented  in  graphical  form  in  the  next 
section.  While  the  process  times  for  these  examples  ranged 
from  one  half  to  seven  minutes,  we  point  out  that  the  DELSIM 
system,  due  to  its  general  nature,  does  not  take  into  account 
the  nice  structure  of  the  system  of  differential  equations  at 
hand.  A program  which  takes  into  account  the  special 
structure  of  the  differential  equations  (3.1.1),  (3.1.2) 
could  handle  much  larger  examples  and  would  also  result  in 
considerable  savings  in  computer  time  and  storage. 

Defining 

5.k(t)  = Gk(t)e , k^O 

we  note  that  the  i-th  entry  of  £k(t)  is  given  by 

gki(t)=P[X(t)=k|X(0)=0,J(0)=i]. 

It  is  easily  seen  from  (2.2.3)  and  (2.2.4)  that 


■ 


(3.3.1  ) 


^ a.0ct )=cn-HCt )T°ApJa0Ct ) 

2.0(0)=^ 


and  for  k>l 


ft  ak(t)  = [T+H(t)T°A0]£k(t)  + {l-H(t)}T0A°3.k_1  (t) 

ak(o)=o. 


\ (3.3.2) 


The  infinite  system  of  differential  equations  above  needs  to 
be  truncated  at  a sufficiently  large  value  of  the  index  k 
before  any  numerical  method  can  be  implemented  to  solve  it. 
It  does  not  appear  tractable  to  develop  any  optimal  methods 
for  such  truncation.  In  the  absence  of  such  criteria  it 
appears  practical  to  truncate  the  system  at  k;Xp+3a  where  a 
is  the  standard  deviation  of  the  stationary  distribution. 
Using  such  methods  we  are  confident  that  the  equations 
(3.3.1),  (3.3.2)  can  be  solved  for  queues  for  which  the 
value  of  Ay  is  even  moderately  large.  But  as  Ay  becomes 
very  large,  say  over  200,  the  systems  of  differential 
equations  above  could  become  stiff  and  pose  considerable 
difficulty  in  solving  them  numerically.  Further  work  along 
these  lines  is  under  way  and  will  be  the  subject  matter  of  a 
forthcoming  paper. 


3.4  Computational  Examples 

To  illustrate  the  computations  of  the  moments  of  the 
queue  length,  we  considered  the  following  five  PH-distribu 
tions  for  the  inter-arrival  times: 


A 


Exponential  (50) , 

(3.4.1) 

E (1 0 , 500), 

(3.4.2) 

0 . 2 E ( 5 , 62.5)  + 0.8E(5,  1000), 

(3.4.3) 

0 . 8E ( 2 , 400 ) +0 . 2 E ( 8 , 100), 

(3.4.4) 

0 . 5E (5  , 1 56 . 25 ) + 0 . 5E (5  , 625)  , 

(3.4.5) 

where  E(n,a)  denotes  the  Erlang  distribution  with  density 

f(x)=iw  e"“x  x"'1-  x=>0- 

While  each  of  these  PH-di stri butl ons  has  the  same  mean  0.02, 

the  distributions  are  qualitatively  very  different  as  the 

graphs  of  their  density  functions  (Figures  1-7)  show.  The 

-4 

variances  of  these  distributions  are  respectively  4x10  , 

0.4xl0“4,  llxlO'4,  10.7x1 0~4  and  2.528xl0'4. 


Figure 
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Figure 
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For  the  service  time  we  considered  the  following  three 


* 


distributions  each  of  which  has  mean  1: 


Exponential 

0), 

(3.4.6) 

R(0,2) 

(3.4.7) 

A discrete  distribution  which  has  mass  0. 
each  of  the  points  0.5  and  1.5. 

5 at 

(3.4.8) 

Figures  8-11 

give  the  graphs  of  the  mean 

(NU1 ) 

and  the 

three  central  moments  ( N U 2 , NU3 , NU4)  respectively  of  the 
queue  length,  plotted  against  time  (X),  for  the  five  PH/G/°° 
queues  each  of  which  has  the  same  service  time  distribution, 
viz..  Exponential  (1),  and  which  have  as  their  respective 
inter-arrival  time  distributions  the  five  PH-di stri butions 
given  by  (3.4.1)  - (3.4.5).  Figures  12-15  present  the 
graphs  of  the  moments  now  under  the  assumption  that  the 
service  time  distribution  is  R(0,2).  Finally,  Figures  16-19 
present  the  graphs  obtained  under  the  assumption  that  the 
service  time  has  the  two-point  distribution  given  by  (3.4.8). 
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An  examination  of  the  graphs  of  the  mean  queue  length 
(NU1 ) given  in  Figures  8,  12  and  16  shows  that  (for  fixed 
service  time  distribution)  each  of  the  five  PH-distributions 
considered  here  results  in  almost  the  same  value  for  the 
mean  queue  length  at  every  point  x.  The  Insensitivity  of 
even  the  time-dependent  mean  queue  length  to  substantial 
random  variability  in  the  arrival  process  may  be  deceptive. 
An  examination  of  the  computed  curves  for  the  time-dependent 
curves  of  the  second,  third  and  fourth  central  moments  show 
that  the  latter  are  all  highly  sensitive  to  variability  in 
the  inter-arrival  times.  In  fact,  increased  variability  in 
the  latter  manifests  itself  in  the  same  qualitative  order  in 
all  the  higher  moments. 

Although,  for  the  initial  conditions  chosen  in  our 
examples,  the  approach  to  "steady-state"  is  very  rapid,  the 
higher  values  of  the  central  moments  can  only  be  explained 
by  a more  erratic  behavior  of  the  path  functions.  These 
observations  indicate  that  sole  reliance  on  simple  analytic 
expressions  for  mean  queue  lengths  (where  these  are 
available)  may  lead  to  serious  errors  in  the  qualitative 
interpretation  of  the  behavior  of  stochastic  models. 


J 


CHAPTER  IV 


THE  PH/PH/®  QUEUE 

4.1  Introduction 

Consider  the  PH/G/®  queue  and  assume  that  the  service 
time  distribution  H(*)  is  also  of  phase  type.  Such  a model 
will  be  denoted  by  PH/PH/®.  In  this  Chapter  we  shall  set  up 
the  differential  equations  for  this  model  indicating  their 
proofs  briefly.  In  the  course  of  our  discussion  we  shall 
also  present  an  alternate  proof  of  an  Interesting  theorem  on 
the  M/PH/®  queue  due  to  M.  F.  Neuts  and  J.  L.  Jain  [8]. 

To  be  specific  let  us  assume  that  the  service  time 

c.d.f.  H ( • ) which  is  of  Phase  Type  has  a representation 
(£,S)  and  consists  of  n phases.  Once  again  we  assume  that 

en+1=0  and  without  loss  of  generality  that  the  representation 

U,S)  is  so  chosen  that  the  PH-Renewal  Process  defined  by  it 
is  irreducible. 

For  the  PH/PH/®  queue  we  are  interested  in  the  random 
variables  X.(t)=the  number  of  customers  in  phase  j of 

J 

service  at  time  t+,  j=l,...,n  and  J(t)=the  phase  of  the 

. 

arrival  process  at  time  t+.  We  let,  for  kA>0,  t>0,  l<i,j<m, 


(4.1.1) 


and  denote  the  m x m matrix  defined  by  these  entries  by 


G(kj  , ,kn  ,t) . For  t>0  and  | z ^ | < 1 , i=l,...,n  we  also  define 
the  generating  function 


G ( z 1 >•••  ,z  ,t) 


k^O,.  . »kn>0 


k k k 

z i ^ z22.  . .znn  G(k1  ....  ,kn,t). 


(4.1.2) 

In  the  sequel  we  shall  simply  write  l£  and  z_  respectively  to 
denote  the  vectors  (kj,...,k  ) and  (z^,...,zn),  and  it  will 
be  implicitly  understood  that  k^O  is  an  integer  for 


i=l,...,n  and  |z. |<1  for  all  i=l,...,n.  We  also  let 


S r ( t ) = { (u1  , . . . ,ur)  : u^O,  zu.<t),  t> 0 , r>l  . 

Further  for  u_eS  (t),  9r(t,u_,k_)  will  denote  the  probability 
of  obtaining  E-j  k-j  times,...,  En  kn  times  in  r multinomial 
trials  each  of  which  can  result  in  any  one  of  the  n+1 
mutually  exclusive  and  collectively  exhaustive  events 
E.|,...,En  + i with  the  trials  having  the  probabilities  given 
by  the  vectors 


(3exp[S(t-Ui  - ... -u.  ) ] , 1 - Bexp[  S ( t-u-j  - ...  - u ..)]•£) ; i = l r 

respectively.  Note  that  3_r  ( t ,u_,k_)  , where  ke<r,  is  the 

conditional  probability,  given  there  are  r arrivals  in  (0,t] 

and  these  occur  at  u^  , u^+u2 , . . . ,u^  + . . .+ur»  that  k^  of  these 

customers  are  in  phase  j of  service  at  time  t+  for  j=l,...,n 
n 

and  the  rest  r-zk.  depart  in  ( 0 , t ] . We  also  have 
1 J 


def  kl  kn 

g (t,u,z)  = i z ' ..z  n g (t  ,u  ,k) 
k>0  1 nr 

ke<r 


r 

= n 

i = l 


{* 


Bexp[S(t-Uj  -...-u^]  • z + 1 -Bexp[S(t 


ur. 


■u.)].ej 


(4.1.3) 


With  these  notations  we  are  now  ready  to  discuss  the  PH/PH/« 
queue . 


4.2  The  Analysis  of  the  PH/PH/°°  Queue 

Lemma  4.2.1:  For  t>0, 

oo  r r 

G(0,t)=exp(Tt)+  E / / n exp(Tu.)*T°A 
r=l  Sr(t)U  = l 1 

exp[T(t-u]  -...-ur)]gr(t  ,u  ,0)du 

and  for  lk>0, 

oo  r r 

G ( Jc , t ) = e / < n 
r=ke  Sr(t)Li=l 

gr(t,u,k)du 

Proof : This  Lemma  is  obtained  by  conditioning  on  r,  the 
number  of  arrivals  in  (0,tj,  and  the  arrival  epochs  u-j  , 
u^  + u2  , • . • ,u-j  + . . . +ur  , and  applying  the  law  of  total 
probability. 

Theorem  4.2.2:  (Basic  System  of  Differential  Equations): 

For  t>*0, 

fj  G(z,t)=[(T+T0A°)+{£exp(St).(z-e)}T0A°jG(z,t)  (4.2.3) 

and 


exp  (Tu  . ) •T°A°^*exp[T(t-u1 -...-ur)  ] 


G(z,0)=I 


(4.2.4) 


Proof : The  above  Theorem  is  obtained  by  using  Lemma  4.2.1, 

(4.1.2)  and  (4.1.3).  We  omit  the  details  which  are  analogous 
to  those  in  the  proof  of  Theorem  2.2.2. 

By  applying  Proposition  (2.2.8)  to  (4.2.3),  (4.2.4)  we 

have 

Theorem  4.2.5:  The  basic  system  of  differential  equations 
for  the  PH/PH/°°  queue  given  by  (4.2.3),  (4.2.4)  has  a unique 
solution  in  [0,°°).  Further  this  solution  is  bounded  in 
[0,o°) . 

Using  Levinson's  Theorem  quoted  as  Proposition  2.4.1, 
we  have 

Theorem  4.2.6:  As  t-*°° , the  solution  G(z_,t)  to  (4.2.3), 
(4.2.4)  converges  to  a matrix  G(zJ  all  whose  rows  are 
identical.  Also  G(e')=0  where  e_'  = ( 1 , . . . , 1 ) eRn  . 

This  implies 

Theorem  4.2.7:  The  process  ^(X-j  (t ) , . . . ,Xn  (t ) ,J  (t ) ) : t>0^ 

is  ergodic  and  has  a stationary  distribution  independent  of 
the  initial  conditions.  This  stationary  distribution  is 
given  by  G [z) . 

Before  we  proceed  with  the  discussion  of  the  moments, 
we  present  the  following  interesting  result  on  the  M/PH/°° 
queue  due  to  M.  F.  Neuts  and  J.  L.  Jain  [8]. 


! 

i 


heorem  4.2.8:  Consider  the  M/PH/®  queue  (which  is  obtained 


by  setting  a=l  and  T=-A  in  the  PH/PH/®  model).  Given 

X-j  ( 0 ) = . . . = Xp  ( 0 ) = 0 , X1  (t) , . . . ,Xn  (t)  are  independent  for  any 

fixed  t>0 , and  for  j = l,...,n,  X,(t)  has  a Poisson  distribu- 

t J 

tion  with  parameter  A/3exp(Sr)e  .dt , where  e.  is  the  n-vector 

0 J J 

whose  j-th  component  is  one  and  all  other  components  are 
zero.  Also  the  process  (X-j  (t) . ,Xn(t) ) is  ergodic,  and  in 
the  steady-state  these  r.v.s.  are  independent  with  X. 

J 

following  a Poisson  distribution  with  parameter  -ApS^e.. 

%J 

Proof : For  the  M/PH/®  queue  we  get  from  (4.2.3)  and  (4.2.4), 

G(z,t)  = ACeexp(St)  • U-e)  }GU,t) 

G(z,0)=l  . 

It  is  easily  verified  that  the  solution  to  the  above  is 
given  by 

n ft 

G(z^,t)=  n exp  -A/Bexp(Sr)* 

j = l L 0 

whence  the  theorem. 

-1 

Remark:  Note  that  -bS  e . = /Bexp(ST)e  .dr  is  the  expected 

-J  0 J 

time  spent  by  a customer  in  phase  j of  service.  In  the 
light  of  this  fact,  the  values  obtained  for  the  parameters 
of  the  Poisson  distributions  in  the  above  theorem  are  very 
intuitive. 

We  now  return  to  the  general  discussion  of  the  PH/PH/® 
queue  and  define  the  factorial  moments 


e.  d-r 

J 


(1 


yi(5)s!E(x!k)(t)|J(0)n1,X1(0)-...-X  (0)=0),  k>l  , l<i<m, 

* \J  sJ 

1 < j<n . 

Also  the  vector 

lLjk)(t)=(yjj)(t),...  ] ( t) ) ' , j = l , . . . ,n , k>l. 

By  differentiating  (4.2.3)  with  respect  to  z.  and  setting 

J 

z=e,  we  can  obtain 


Theorem  4.2.9:  For  t>0,  1 < j < n , 


^ )(t)=(T+T°A°)1Lj1  ^ ( t )+{  6.eSte.j  }T° 


41  )(°)  = 0 
J 


(4.2.10) 


and  for  k>2 


^ y.jk^(t)  = (T+T0A°)yjK;(t)  + k{£eJ  Le^  }T°A%  j 


(k-1) 


yjk)(0)=0. 


(4.2.11  ) 


Further  as  t-*°°,  for  any  k > 1 , l<j<n,  v}.  '(t)  converges  to  a 

J 

finite  vector  all  whose  components  are  equal. 

We  conclude  our  discussion  by  presenting  the  following 


Corollary  4.2.10:  For  1 < j <n , as  t-»-°°,  y:  ;(t)  converges  to  a 

J 

vector  each  of  whose  components  is  equal  to  -AjbS  ^e. , where 

J 

A-1  is  the  mean  inter-arrival  time  given  by  A'^-aT^e. 
Remark : Noting  that  -£S_1e.  is  the  expected  time  spent  by  a 

J 

customer  in  phase  j of  service,  we  see  that  the  above  result 
is  very  intuitive. 


Proof : From  (4.2.10), 

ePiP)(t)  = 0l°/(£eSTe.)dT 
- 3 0 J 

-*■£1°  ( "Is -1  e . ) as  t-*».. 

J 

=-xes" 1 , 

for,  £=-XaT_1  and  T°=-Te. 

The  proof  is  complete  by  appealing  to  the  last  statement  of 
Theorem  4.2.9. 
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